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Abstract. We give a characterization for the extreme points of the convex set of correlation 
matrices with a countable index set. A Hermitian matrix is called a correlation matrix if it is 
positive semidefinite with unit diagonal entries. Using the characterization we show that there 
exist extreme points of any rank. 
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1. Introduction 

Let X be a fixed index set, and let F denote the field of either real or complex numbers. A 
mapping C : X x X — >• F is an (auto) correlation function if it is positive semidefinite, Hermitian 
(or symmetric when F = R), and C(i, i) = 1 for all % G X. This term comes from the theory of 
stochastic processes. It the case where X is either finite or countable, such a function is called a 
correlation matrix. The set of correlation functions is convex. The problem of determining the 
extreme points of this set has been studied extensively, but the research is mostly concentrated 
on finite matrices [1, H, E] • 

Christensen and Vesterstr0m [2] considered the complex case and proved that, when the 
cardinality of X is greater than 3, there exists a rank 2 extreme point. Moreover, they showed 
that when X is finite with n elements, the rank of any extreme point is at most y/n. Loewy [9| 
proved that, in fact, for any r < \/n there exists a rank r extreme point. Grone, Pierce, and 
Watkins [4] showed that, in the real case, a similar result holds but the criterion for the rank 
is then r 2 + r < 2n. Finally, Li and Tarn were able to give a simple characterization for the 
extreme points in the case where X is finite [HI Theorem 1 (b)]. This result holds in both real 
and complex cases. 

The complex correlation matrices are used also in quantum mechanics where they appear, 
e.g., as the structure matrices of certain T-covariant observables 0, CD, E]: Finite correlation 
matrices occur in the context of angle observables while the case of countable X is associated 
with covariant phase observables (when X = {0, 1, 2, ...}) and with box localization observables 
(when X = Z). In these cases, the extreme correlation matrices correspond to the extreme 
T-covariant measurements; accordingly, Holevo stressed the importance of characterizing them 
[5]. For finite X, the solution is provided by Li and Tam as mentioned above. Their result has 
been generalized to other finite-dimensional covariance systems by D'Ariano [3j. 

The purpose of this note is to give a characterization of extreme correlation matrices in 
the case of countable index set (which can be chosen to be N without restricting generality). 
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We consider both real and complex cases. The result is given by Theorem 1 below. It is a 
generalization of Theorem 1 (b) of [8j. We also show that there exist exteme matrices of any 
rank (G N U {oo}). 



2. Notations and basic definitions 

Let M. be the set of infinite Hermitian (i.e. conjugate symmetric) F-valued matrices indexed 
by N = {1,2,3, ...}, considered as a real linear space with respect to the usual entrywise sum 
and scalar multiplication. Let V be the vector space over F consisting of F-valued sequences 
c = (c n ) nS N such that c n ^ only for finitely many n G N. For any fceN, let G V be the 
sequence (5kn)nm where 5k n is the Kronecker delta (so that (ekjkeN forms an algebraic basis 
of V). Each matrix M = (M nm ) n me N G M. defines a Hermitian sesquilinear (or symmetric 
bilinear when F = R) form M : V x V — > F via 

M(c, d) = } j c^M nm d m , c,deV. 

n,m£N 

(Notice that the above sums are finite.) Let C C M. be the convex set of Hermitian positive 
semidefinite matrices with unit diagonal elements, i.e. the set of correlation matrices. Recall 
that C G M. is positive semidefinite if and only if C(c, c) > for all c G V. 

Let £ X (N) be the Banach space (over F) consisting of F-valued sequences c = (c„) nG N such 
that EneNM < 00 ■ We denote the ^-norm 2^„ eN |c„| of c G £ X (N) by \\c\\ x . Let £ X (N) X 
denote the topological antidual of ^(N), i.e. the linear space (over F) of continuous antilinear 
functionals £ X (N) — > F. We equip £ X (N) X with the standard operator norm 

/->sup{|/(c)| \ce£\N), ||c||a<l}. 

Naturally, £ 1 (N) X is isomorphic with £°°(N) = {d : N -> F | sup nG pj \d n \ < oo} but we do not 
actually need this fact. 

Let C G C. Since C nn — 1, n G N, and the principal 2 x 2-minors of C are nonnegative, we 
have sup nmgN \C nm \ < 1, implying that 



sup 



^2c nm c m < \\c\lx, cef(N). 



meN 



Hence, C defines a continuous linear map C : £ 1 (N) — * £ 1 {N) X by 

[C(c)] (d) = d~nC nm c m , c,de £\N). 

n,mSN 

where the double series converges absolutely. 
The rank of C G C is defined by 

rankC = dimC'(£ 1 (N)) G N U {oo}, 

i.e. rankC is the dimension of the linear space C(£ 1 (N)) C £ 1 (N) X . 

3. Technical lemmas 

Before we can prove the main result of this note (Theorem [1]) , we need some information on 
the structure of a matrix C G C. 

The matrices C G C have the following characterization. Let C G M.. Then C G C if and only 
if there exists a separable Hilbert space TC over F and a sequence of unit vectors (r] n ) neN C TC 
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such that C nm = (rj n \rj m ) for all n, m G N [lj (see also [10], Exercise 8.7]); here ( • | • ) denotes the 
inner product of Ti. (linear in the second argument) and later we use the symbol || • || for the 
norm of H. In the above references, this result is given only in the complex case, but the real 
case follows easily by noting that a complex Hilbert space (TC, ( ■ \ ■ )) can be interpreted as the 
real Hilbert space equipped with the inner product Re ( • | • ). Note that TC and the sequence 
(jin) are not uniquely determined by C . 

Let C G C, and choose a sequence (r) n ) of unit vectors in a Hilbert space TC (over F) such 
that C nm = {rj n \rj m ) for all n,m G N. These choices will remain fixed in this section. For each 
c G ^(N), the series J2neN c nVn converges (absolutely) in TC, so we can define a linear map 
$:^(N)->Wby 

(1) $(c) = ^c n r/ n . 

neN 

Since ||$(c)|| < ||c||i the map $ is continuous. 

Next define another linear map $* : TC — > £ 1 (N) X by 

[$*(y)](c) = X;^nb). ^W, ce£\N), 

neN 

where the series converges absolutely. The map $* is continuous since 

sup {|[$»](c)| I c G £\N), Helix < 1} < |M|. 
It is clear from the above formulas that we have 

(2) [*»](c) = <*(*>, pen, cel\N). 

Lemma 1. (a) lin F {r? n | n G N} = $(^(N)) = (ker$*) ± . 

(b) C = 

(c) C(^(N)) = 

(d) rankC = dim$*(ft) = dim^^N)) = dimlin F {r/„ | n G N}. 

Proof. Since i] k = $(e fc ) by JT]), we have lin F {^ n | n G N} C $(^(N)). On the other hand, 
the series in (JTJ) converges in norm, so $(£ 1 (N)) C lin F {^ n | n G N}. This proves the first 
equality in (a). It follows from p that $(^(N)) C (ker^*)^ and hence $(^(N)) C (kerS*) 1 - 
(the orthogonal complement is closed). Assume now that (p G (ker$*)^ fl $(£ 1 (N)) ± . Then 
[$*(</?o)] (c) = ($(c)\(p ) = for all c G by ©, so y? e ker$* and hence y? = 0. Now (a) 

is proved. 

By dH), the continuity of $*, and (j2J), we get 

[**(*(c))](<0 = ^c n [$*(r /n )](rf) = ^c n ($(rf)|r ?n ) 

n6N neN 

= ^Y,^ m \Vn)Cn=[C(cj\(d), C,de£\N). 

neN m€N 

This proves (b). 

To prove (c), we first note that C (l* (N)) C $*(H) by (b), so C(^(N)) C <£>*(H). On the other 
hand, by using the relation $(£ X (N)) = (ker^*)^ (see (a)) and the fact that $* is continuous, 
one gets 

= $*($(£!(N))) C $*($(^(N))) = C(f-(N)), 
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and hence $*(7Y) C C^N)), and (c) is proved. 

It remains to prove (d). By (a), we have linF{?7 n | n G N} = $(€ 1 (N)) = /C, where /C = 
(ker$*) ± . Since any finite dimensional space is closed, \m^{r] n n e N} = dim$(£ 1 (N)) = 
dim/C. Now $*k(/C) = since $*(/C x ) = 0. In addition, the map : K -> is 

a linear bijection, so dim/C = dim$*(7Y). It follows from (c) that dim(7(£ 1 (N)) = dim$*(7i), 
i.e. rankC = dim$*(7i). This completes the proof of (d). □ 

Remark 1. From Lemma 1, one sees that, without restricting generality, we may assume that 
dim Ti = rankC and the sequence (r?n)neN spans Ti. In particular, if rank C = r < oo then we 
may assume that Ti = F r . 

For each cp G Ti, we let \^){^\ denote the bounded operator ip i— > on 7i. We use the 

symbols £ s (Ti) and T S (TC) for the real Banach spaces of selfadjoint bounded and selfadjoint 
trace class operators on Ti, respectively. The following lemma is well known in the case of a 
complex Hilbert space. However, real Hilbert spaces are not so frequently used in the literature, 
so we give the proof here for the real case. 

Lemma 2. The topological dual ofT s (Ti) is isomorphic to C s (Ti). 

Proof. Assume that Ti is real (i.e. F = R). For each R G C s (Ti), define a linear mapping 
F R : %{Ti) -> R by F R (T) = tr[RT}. The functional Fr is continuous, since \F R (T)\ < 
\\R\\\\T\\ tT where || • || and || ■ || tr are the operator and trace class norms, respectively [3, 
Theorem 4.1.4 (3)]. Let then F : %{Ti) — > R be a continuous functional. Define an R-linear 
map F on the whole trace class T(Ti) of Ti by F(T) = \F(T + T*), T G T(H). Since 
||T*|| tr = ||T|| tr , it follows that |F(T)|| < ||F||||T|| t r and F is continuous. Hence, there exists a 
bounded operator S such that F(T) = tr[ST], T G T(Ti), and ||F|| = \\S\\ Theorem 4.1.4 
(3)]. For any T G %(Ti) we have F(T) = F(T) = |tr[(5 + S*)T\. (Note that tr[A*} = ti[A], 
A <ET (Ti).) Define R = \{S + S*) G C s {Ti) so that F = F R . The mapping R i-> F fl is isometry. 
Indeed, = ||i?|| since tr[i?|(T + T*)] = tr(i?T), T G T(W), and |||T + T*|| tr < 1 when 

imi tr <i. ' ' □ 

4. Extreme points of C 

The following theorem characterizes the extreme points of the convex set C of F-valued 
correlation matrices. The characterization is an infinite-dimensional analogue of [H Theorem 1 
(b)]- 

Theorem 1. Let C G C, and let (r] n ) be a sequence of unit vectors in a Hilbert space Ti (over 
¥) such that C nm = (r] n \r] m ) for all n, m G N. Denote JC = \m ¥ {r] n \ n G N} and 

S = lin M { | r] n )( Vn | | n G N} C %{K). 

Then C is an extreme point of C if and only if S is trace norm dense in %(JC). 

Proof. Define $ and $* using (j} n ) as in the preceding section so that Lemma Q] holds for them. 
Notice that K = (ker$*) ± = $(^(N)) CHby Lemma U (a). 

Assume that S is not dense in %{K). Since the dual of T S (JC) is C S (JC) by Lemma 2, it follows 
from the Hahn-Banach theorem that there is an R G C S (K) such that R ^ and tr[i?T] = for 
all T eS. Define B : i\N) -> £ 1 (N) X by B = &*\ K R$. Now B ^ 0, because $*|^ is injection 
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and $(£ 1 (N)) is dense in fC by Lemma [T] (a) . Let B nm = [J3(e m )l (e n ) for all n,m G N. Then 
by © and Q, we get 

B nm = [$*{m{e m ))]{e n ) = ($(e n )| J R$(e m )) = (rj n \Rr] m }, n,meN. 

Since i? is selfadjoint, it follows that the matrix (B nm ) is Hermitian. In addition, B nn = 
tr[R\rj n ) {r] n \} = for all n G N, since \r] n )(r] n \ G S for all n G N. Let e = ||-R|| _1 (where ||i?|| is 
the operator norm of R). Since R is a selfadjoint bounded operator, I ± eR > where J is the 
identity operator of /C. Hence, for all c G ^ 1 (N), one gets from ([2]) that 

(C7±eS)(c,c) = [($*$ ±e$*i?$)(c)](c) = ($(c)|(I± efl)$(c)) > 0. 

It follows that the matrix C ± eS is positive semidefinite. Since, in addition, B nn = for all n, 
we get C ± eB G C. Clearly C = \{C — eB) + |(C + e5), so C is not an extreme point. 

Assume then that C is not an extreme point. Then C = \{C\ + C 2 ) for some C x , C 2 G C, 
where C x ± C 2 . Define £ = \{C X - C 2 ). Then B ^ 0, and B e M, with = for all n G N. 
Since C ± B G C, one gets C(c, c) ± _B(c, c) > for all c G V, and hence 

(3) \B(c,c)\ <C(c,c), cGV. 
Next we show that 

(4) if c G ker $ n V and d G V then E(c, d) = 0. 

To that end, let c G V be such that $(c) = 0, and let d G V. Then C(c, c) = [$*$(c)j (c) = 0, 
so also (C + £>)(c, c) = by ((3|). Since C and C + B are positive semidefinite and Hermitian, 
the Cauchy-Schwarz inequality gives C(c, c?) = and (C + S)(c, d) = 0. Hence also B(c, d) = 0, 
proving (TjJ. 

It follows from (TjJ that B(c x ,d x ) = B(c 2 ,d 2 ) whenever c x , c 2 , d x , d 2 G V are such that c x — 
c x ,d x — d 2 G ker$. (Notice that B is Hermitian.) Hence, we have a well-defined sesquilinear 
form R : $(V) x $(V) — ► F defined by R(ip, ip) = B(c, d) where <3>(c) = and = ?/>. Since 
B is a Hermitian, also R is such. 

Let now ip G 3>(V), and let c G V be such that p = $(c). Then by ([3|), Lemma CD (b) and 
0, we get 

\R(<p,<p)\ = \B(c,c)\ < C(c,c) = [**(*(c))](c) = <$(c)|$(c)> = |M| 2 , 

so the polarization identity gives 

svp{\R(<p,4>)\ \<p,?P e $(V), |M| < 1, HVII < 1} < 00. 

Since 3>(V) = lin F {^ n | n G N} is dense in K, by Lemma Q] (a) it follows that there is a bounded 
operator R : /C — > K., such that R(p,ip) = (p\Rip) for all p,ip E $(V). Since i? is Hermitian, 
-R G £ S (/C). Now $(e„) = r] n for all n G N, so 

tr[ J R|r]„)(?7 n |] = (r} n \Rr} n ) = B(e n , e n ) = B nn = 0, n G N. 

This implies that tr[_RT] = for all T G S. 

Now if S is dense in T 8 (K), it follows that ti[RT] = for all T G T S (.M), implying .R = 0. 
But then 

B nm = B(e n , e m ) = R(r} n , r] m ) = (r} n \Rr] m ) = 0, n, m G N, 
which is impossible, since B 7^ 0. Hence, 5 is not dense in %{K,). The proof is complete. □ 

The following proposition shows that there exist extreme points of C of any rank (g NU{oo}). 

5 



Proposition 1. Let Af be a (separable) Hilbert space over F, and let {<p n \ n G iV} CJV \ {0} 
be a countable set dense in Af ' . Define 



Then the matrix C 



C nm =\\f n \\ 1 \\Vm\\ 1 ( ( Pn\ l Pm), n, m G N. 

(C nm ) is an extreme point of C with rankC = dim Af. 



Proof. Let r\ n = \\(p n \\ Vn for all neN. Since the set {<p n n G N} is dense in Af, it is 
clear that lmw{r] n n 6 N} = Af. Therefore, rankC = dim Af by Lemma [Q (d). We proceed to 
show that the set S = \m K {\r] n ) (r] n \ \ n G N} is trace norm dense in T s (Af). To that end, let 
T G %{Af), T^O, and e > 0. Using the spectral representation [7], we get 

t = y^t n |0 n )(0 n |, 

neN 

where t„6R for all fiGN, and (0 n )neN is an orthonormal sequence in Af, the series converging 
in the trace norm II ■ |L. Choose n G N such that 



(5) 



T 



no 

E 

n=l 



tn^n) (<t>n\ 



Then, for each n = 0, ...,n , pick fc n eN such that \\(f> n — fk n \\ < e ( 6 1 1 CZ -1 1 1 tr ) 1 - This is possible 
because {y? n | n G N} is dense in Af. Now we use the fact that for ifi,ip2 G A/", such that 
ll^ill = 1 and 11^1 — ^2 II < 1, we have 



Wl)tyl\- 1^2)^2 



<3| 



^2 



This can be proved [HI the proof of Lemma 5] e.g. by using the duality C(Af)* = T(Af) [3, p. 
60] where C(Af) denotes the set of compact operators. Applying this result, we get 



no 



(6) 

Now (p2) and 

(7) 



^2t n \(f) n )(<j) n \ -^2t n \tp kn )(ip kn \ ^ < \ f 



n=l 



n 



n=l 



n=l 



imply that 



»o 



n=l 



< e. 



Since ip kn = \\<fk n \\Vk n we have YTn^nWkf) (<£k n \ G 5. Thus 5 is dense in T s (Af). Now C is an 
extreme point by Theorem 1. The proof is complete. □ 

The following remark demonstrates that there exist extreme correlation matrices which can- 
not be constructed by using Proposition 1. 

Remark 2. It is easy to generalize the construction of finite rank extreme points of Li and 
Tarn P, Section 2.1] to the infinite dimensional case. 

First we consider the complex case. Let r G N and choose Ti = C r . Let (f n ) n= i be an 
orthonormal basis of C r . Then the basis of the finite dimensional space T s (C r ) = £ s (C r ) consists 
of operators \ f n ) (/ n |, |/ n ) (f m \ + \ f m ) (f n \ and i \f n ) (f m \ -i\f m ) (f n \ where n,m G {1,2, ...,r} 
and n < m. Define (r^-sequence as follows. Enumerate the unit vectors /„, 2 _1//2 (/„ + f m ), 
2- 1 / 2 (f n 4- n,m G {1, ...,r}, n < m, to get the first r 2 vectors 7] n and define = f\ for 
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all n > r 2 . Thus we have a sequence (^ n )ngN of unit vectors and we can define a C e C by 
C nm — (rj n \rj m ). Since, for example, 

i l/n) - 1 (/n| = |/n) (fn\ + \fm) (fm\ ~ 2 • |2- 1/2 (/„ + i/ m )> ^"V 2 ^ + i/ m ) | 

it is easy to see that \in^{\ri n ) {r] n \ \ n G N} = T s (C r ). Hence, C is extreme. Moreover, 
rankC = dimlinc{r/ n n 6 N} = r. 

Let {fin) be any sequence of unit vectors (in some Hilbert space) such that C nm = (fj n \fj m ). 
Since (fj n \fj m ) = 1 for all n,m > r 2 it follows that f) n = fj r 2 +l for all n > r 2 by Cauchy-Schwarz. 
Hence, there exist only finitely many different vectors fj n , so the sequence (fj n ) cannot be dense 
in the surface of the unit ball. It follows that the sequence (fj n ) cannot be obtained by the 
method of Proposition 1. 

In the real case we take 7i = W with an orthonormal basis (f n ) n= i- Now the operators 
\f n ) (f n \ and \f n ) (f m \ + \ f m ) (/ n |, n,m G {1,2, ...,r}, n < m, generate T S (W) and we can 
choose the first r(r + l)/2 ^-vectors to be /„ and 2~ 1 ^ 2 (f n + f m ), n < m. Then we set r/ n = fi 
for all n > r(r + l)/2. As before, we notice that the matrix C nm = (rj n \rj m ) is an extreme point 
of rank r and it cannot be obtained by Proposition 1. 

Hence, to construct a finite rank extreme, Proposition 1 gives unnecessarily many different 
vectors rj n . 
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